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The 15 N(p,7) 16 reaction provides a path from the CN cycle to the CNO bi-cycle and CNO 
tri-cycle. The measured astrophysical factor for this reaction is dominated by resonant capture 
through two strong J w = 1~ resonances at En — 312 and 962 keV and direct capture to the ground 
state. Recently, a new measurement of the astrophysical factor for the 15 N(p, 7) ie O reaction has 
been published [P. J. LeBlanc et ai, Phys. Rev. C 82, 055804 (2010)]. The analysis has been 
done using the ii-matrix approach with unconstrained variation of all parameters including the 
asymptotic normalization coefficient (ANC). The best fit has been obtained for the square of the 
ANC C 2 = 539.2 fm -1 , which exceeds the previously measured value by a factor of w 3. Here 
we present a new i?-matrix analysis of the Notre Dame-LUNA data with the fixed within the 
experimental uncertainties square of the ANC C 2 = 200.34 fm -1 . Rather than varying the ANC we 
add the contribution from a background resonance that effectively takes into account contributions 
from higher levels. Altogether we present 8 fits, five unconstrained and three constrained. In all the 
fits the ANC is fixed at the previously determined experimental value C 2 = 200.34 fm -1 . For the 
unconstrained fit with the boundary condition B c — S C (E2), where E2 is the energy of the second 
level, we get S(0) = 39.0 ± 1.1 keVb and normalized x 2 = 1.84, i.e. the result which is similar 
to [P. J. LeBlanc et ai, Phys. Rev. C 82, 055804 (2010)]. From all our fits we get the range 
33.1 < 5(0) < 40.1 keVb which overlaps with the result of [P. J. LeBlanc et al, Phys. Rev. C 82, 
055804 (2010)]. We address also physical interpretation of the fitting parameters. 

PACS numbers: 26.20.-f, 25.60.Tv, 24.30.-v, 29.85.-c 



I. INTRODUCTION 

The 15 N(p, 7) 16 reaction provides the path to form 16 in stellar hydrogen burning, thus transforming the CN 
cycle into the CNO bi-cycle and CNO tri-cycle. In stellar environments, the 15 N(p, 7) 16 reaction proceeds at very 
low energies, where it is dominated by resonant capture to the ground state through the first two interfering J™ = 1~ 
s wave resonances at En = 312 and 964 keV, where Er is the resonance energy in the center of mass (c. m.) system. 
There is also a small contribution from direct capture to the ground state of 16 O, which turns out to play an important 
role due to the interference with the resonant amplitudes. 

In our previous paper pQ, the measurement of the asymptotic normalization coefficient (ANC) for the 16 — > 15 N+p 
has been reported. This ANC has been used to fix the non-resonant contribution to the 15 N(p, 7) 16 capture and 
we found that it was impossible to fit the low-energy data from [2]. Moreover, we underscored that to fit these 
experimental data one needs to increase the ANC almost by on order of magnitude. Our calculated astrophysical 
factor using the two-level, two-channel i?-matrix approach led to 5(0) = 36.0 ±6.0 keVb, which is significantly smaller 
than 5(0) — 64 ± 6 keVb reported in [2] but in agreement with the older measurements in [3]. Correspondingly, we 
have found that for every 2200 ± 300 cycles of the main CN cycle, one CN catalyst is lost due to the 15 N(p, 7) 16 
reaction, rather than 1200 ± 100 cycles determined from data of Ref. [2J. Our results were confirmed later by 
[3]. Recently, two new measurements of the astrophysical factor for the 15 N(p, 7) 16 have been published. The 
first measurement was performed at the LUNA underground accelerator facility at the Gran Sasso laboratory [5]. 
This measurement covered only the low-energy region, E < 230 keV, where E is the relative 15 N — p energy. The 
second study of 15 N(p, 7) 16 0, which has been just reported in [B], was performed over a wide energy range at the 
Notre Dame Nuclear Science Laboratory (NSL) and the LUNA II facility. The obtained 5(0) = 39.6 ± 2.6 keVb 
is in a perfect agreement with our prediction 5(0) = 36.0 ± 6.0 keVb pQ. However, in the _R-matrix fitting of the 
experimental data the ANC was used as an unconstrained fitting parameter and the best fit with the normalized 
X 2 = 1-80 has been achieved for the square of the ANC C 2 = 539 ± 138 fm -1 , which is significantly higher then our 
measured value C 2 = 192 ± 26 fm" 1 . The ANC is a fundamental nuclear characteristics [7] and, if it is available from 
independent experimental measurements, it can be varied only within experimental boundaries. An unconstrained 
variation of the ANC to achieve the best fit might signal that some physical input is missing in the reaction model 
and this incomplete knowledge is compensated for by adopting an unphysical value of the ANC. In the case under 
consideration, in the i?-matrix approach the ANC determines the overall normalization of the non-resonant radiative 
capture amplitude and the channel (external) part of the radiative width amplitudes of the both resonances involved. 



Although these amplitudes are small and any sizable impact on the astrophysical factor can be achieved only by a 
significant variation of the ANC, the contribution of the non-resonant amplitude increases toward low energies, which 
is the region of the astrophysical interest. However, not every best fit can be accepted if physical parameters of the 
fit exceed the previously well established experimental limits. A similar situation has occurred in the analysis of the 
14 N(p, 7) 15 reaction for transition to the ground state, where the non-resonant capture contribution is also controlled 
by the ANC for the 15 0(gr.st.) -> 14 N + p Ej. The best fit is achieved at C m 11 fin" 1 / 2 , but this fit was not 
accepted because the recommended ANC value is 7.4 ± 0.5 fm~ 1//2 [8, 9 . The recommended in [5] 5(0) = 0.27 ± 0.05 
keVb factor corresponds to the fit with the recommended ANC but the expanded uncertainty is not due to the best fit 
at higher ANC but to a possible contribution of the capture to the channel spin / = 1/2, which interferes with the 259 
keV resonance. It is an instructive example when missing physics is taken into account rather than an unconstrained 
variation of the ANC, which demonstrates that it is important to keep track on the boundaries of variation of the 
fitting parameters in accordance with the previously available information about them obtained from other sources. 
For example, in the case under consideration these additional parameters are the observable partial resonance widths 
in the proton and a-channels, which have been determined from the i?-matrix analysis [4j [10| of the 15 N(p, a) 16 C 
data [TT] . It is especially important because eventually the results of the fit are aimed to deliver vital astrophysical 
information and to make it reliable all available nuclear physics information should be invoked. 

In this paper, we present our own fits of the data from [6] with the fixed experimental ANC but adding the 
contribution from a background resonance. Results of our fit are compared with the one presented in [6] obtained 
using the AZURE iJ-matrix code [IS]. Besides, we compare the parameters of both fits and present also the fit where 
observable particle widths are consistent with the previous fits of the 15 N(p, a) 16 C reaction. 



II. ANC 



We start with the physical meaning of the ANC. The residue of the scattering 5 matrix in the corresponding 
bound-state pole is expressed in terms of the ANC [13], what provides a model-independent definition of the ANC. 
In the i?-matrix approach the ANC determines the normalization of the external non-resonant radiative capture 
amplitude and the channel radiative reduced width amplitude 14]. In a single-particle approach, the nucleon ANC 
can be expressed in terms of the spectroscopic factor and the single-particle bound-state wave function of the nucleon 
calculated in some adopted mean-field: 

C 2 = 56 2 , (1) 

where 5 is the spectroscopic factor and b is the single-particle ANC, i.e. the amplitude of the single-particle bound- 
state wave function. Note that the isospin Clebsch-Gordan coefficient is absorbed in the spectroscopic factor. In 
such an approach we can consider the ground state of 16 as the bound state ( 15 Np) with the proton occupying the 
single-particle orbital lpi/2- The spectroscopic factor shows the probability of this configuration in the ground state 
of 16 0. Due to the identity of nucleons this probability can be larger than one. In a simple independent particle shell 
model the spectroscopic factor of the lpi/2 state is equal to the number of protons occupying this orbital, i.e. 2. To 
determine the spectroscopic factor from Eq. (TJ) one needs to determine the proton bound-state wave function. To 
this end we adopt the Woods-Saxon potential with the standard geometry, tq = 1.25 fm and diffuseness a = 0.65 
fm and the spin-orbit potential depth of 6.39 MeV. Assuming that the proton is in the P1/2 orbital we obtain the 
single-particle proton ANC b sa 9.96 fm -1 / 2 . If one adopts the ANC from [6], we obtain the spectroscopic factor 
5 = 5.44 versus 5 = 2.1 obtained from our C 2 — 192 fm -1 . Even if we adopt an unrealistically large radius r = 1.50 
fm and a = 0.65 fm for the Woods-Saxon potential, we obtain b = 13.6 fm -1 ' 2 and too high spectroscopic factor 
5 = 2.81. 

Definitely such a high spectroscopic factor obtained from the ANC adopted in |6j requires a physical interpretation. 
In this aspect it would be useful to present the phenomenological spectroscopic factors obtained from the analysis 
of different reactions. For example, the spectroscopic factor deduced from the analysis of 16 0(e, e' p) 15 N reaction is 
5 = 1.27±0.13 |T21[T2]. The proton bound-state wave function deduced from the (e, e' p) reaction is reproduced by the 
Woods-Saxon potential with the geometry ro = 1.37 fm and a = 0.65 fm. The single-particle ANC of this bound-state 
wave function is b — 11.62 fm" 1 / 2 . Using the upper limit of the deduced spectroscopic factor 5 = 1.40 we obtain 
the square of the proton ANC in the ground state of 16 C 2 = 189 fm -1 , which is in a perfect agreement with our 
result pQ. In our paper [1], we presented the spectroscopic factors extracted from different 15 N( 3 He, d) 16 reactions 
including our result. Besides our spectroscopic factor 5 = 2.1 [JJ, three other measurements (see references in [JJ) 
gave 5 < 1.76. The DWBA reanalysis of the 16 0( 3 <i, He) 15 N reaction [TT] performed in [T5] using the proton bound- 
state wave function obtained from the 16 0(e, e'p) 15 N reaction with tq = 1.37 fm and a = 0.65 fm gave even lower 
spectroscopic factor 5 = 1.02. With this low spectroscopic factor we get C 2 ~ 140 fm -1 . Since the measurements in 
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[T7] were done in 1967, the accuracy of the absolute normalization of the differential cross section in [T7] might be 
questionable. Spectroscopic factors below 2 have also been obtained in microscopic calculations [HI [19]. Concluding 
this discussion we cannot find any justification for such a high value of the ANC used in [5J because it leads to an 
unphysical spectroscopic factor. We only can assume that a broad variation of the ANC beyond of the experimental 
limits has been used in [5J to compensate for missing mechanisms in the reaction model. 



III. R MATRIX 



To underscore the role of the ANC in the fitting experimental data in this section we present the expression for 
the astrophysical factor in the R- matrix approach, which we use for the analysis of the experimental data [6J. It is 
two-level, two-channel i?-matrix, which includes the coherent contribution from two 1 _ resonances and non-resonant 
term describing direct capture. The ANC determines the normalization of the non-resonant capture amplitude, which 
describes the external direct capture in the i?-matrix approach and the channel radiative width amplitude, which 
are important for the fitting. But, in addition to [BJ, we add the coherent contribution from a background resonance 
rather than varying the ANC. We can assume that this background pole takes effectively into account contributions 
from distant levels with J 71 " = 1~. The expression for the astrophysical factor in the i?-matrix method for the case 
under consideration can be written as [201 EU 



SW(keVb) = ^^?±^m.#*«viO[ £ (Tl^ nt) ±Tl^ xt )[A-% T r^ ± Mnc+M BG ], (2) 
2 J X J A m x m A ^ > 



V,T=1,2 



where Ajv = 0.2118 fm is the nucleon Compton wave length, 931.5 is the atomic mass unit in MeV, Zj and rrij are the 
charge and mass of particle j and /x^ is the reduced mass of particles i and j, r/ is the Coulomb parameter in the initial 
state of the reaction, Jj is the spin of particle j and Jr is the spin of the resonance, J = 2 ,/+ 1, k 1 = (E + e)/(hc) is 
the momentum of the emitted photon expressed in fm -1 , E is the relative p — A energy, A is the target, e is the proton 
binding energy of the bound state (Ap), L is the multipolarity of the electromagnetic transition, A is the standard 

level matrix for the two-channel, two-level case [20] . = F\ c (k c , r c ) j TC , Jrc is the reduced width amplitude 
for the level r = 1, 2 in the channel c = p,a, P; c (fc c , r c ) is the barrier penetrability factor in the channel c, l c is the 
orbital angular momentum of the resonance and k c is the relative momentum of the particles in the channel c, r c is 
the i?-matrix channel radius in the channel c, rj/^int) = \/2 k^ +1 ^ 2 j v 7 (i nt ), 7j, 7 (mt) is th e internal radiative reduced 
width amplitude for transition from the level v to a bound state (the ground state in the case under consideration), 

7f 7cxt , 7^ Text is the complex channel (external) radiative reduced width amplitude for transition 
from the level v to a bound state given by the expression 



„ lie 2 , 931.5 r , 1/2 T r Z v , sT Z A , (L + 1)L 1 1/2 / " r 

= c V 2 hc Xn ^t r * +/ ^4 + 4 ] V —r- ^ r -WW 

[Fl(k p ,r p ) + G? p (*fc,r p )]W 1/+ i /2 (r p ) < l p L 0\l f > U(Ll s JJ- l p J f ) 



:[/nti + i ^l^fe^ j^C^/nfa], (3) 



F 2 (k P ,r p ) f r L W- n „ lf+ i /2 (r) F lp (k p ,r) 



Int 1 = F { ^ rp) P lp {k p r p ) f dr 



r^ +1 W_ Vftlf+1/2 (r p ) F lp (k p ,r p ) 



^ (vr 'W) ^ r-^:>w v , (4) 



k P r p p P' p J r L+i W_ 7)f j f+1/2 (r p ) Gi p (kp,r P Y 



, r L W_ Vf , lf+1/2 (r) F lp (k p ,r) 7 r L W^^ lf+l/2 {r) G lp (k p ,r) 
mz 2 - / ar^—r— -— - ; ar—r—r— -—— [t>) 



r^ +1 W_ nfjlf+1/2 (r p ) Fi p (k p ,r p ) J r ^ +1 W- Vfilf+1/2 (r p ) Gi p (k p ,r P Y 

Fi p (k p ,r) and Gi p (k p ,r p ) are the Coulomb regular and singular solutions, W- rifi i f +i/ 2 (r) is the Whittaker function 
describing the radial dependence of the tail of the bound state wave function to which transition occurs after the photon 
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is emitted, r/f is the Coulomb parameter of the bound state and If its orbital angular momentum; < l p L 0| //0 > 
is the Clebsch-Gordan coefficient and U(LlfJiI; l p Jf) is the normalized Racah coefficient, Ji is the total angular 
momentum of the system p + A in the initial state of the radiative capture process, which is equal to the resonance 
spin, Ji — Jr, / is the channel spin. The non-resonant capture amplitude describing the external direct capture in 
the i?-matrix method is given by 



M DC = 2C 



e \J n - 5 (k r) L +^u L \ Zp +( 1)L ZA h (L + l)L 1 / 

W^ fJf+1/2 (r p )F(k p ,r p )G(k p ,r p ) < l p L 0\l f > U{Ll f JJ- l p J f )Int 2 . (6) 

As we can see the channel radiative width amplitudes and the direct capture amplitude are proportional to the same 
ANC because both describe the peripheral processes contributed by the tail of the overlap function whose amplitude 
is the ANC. Besides, the channel radiative width amplitude 7„ 7 ( cx t) contains the proton reduced width amplitude 
7,,p. Hence the relative sign of 7„ 7 ( GX t) and Mjjc depends on the sign of 7 I/p . It is also worth mentioning that 7„ 7 (i n t) 
is a fitting parameter and 7„ 7 ( e xt) is a complex quantity because it contains an imaginary part. The radiative width 
is given by equation 

r ^7 = l r / 7 (int) _ r 7/7(oxt)| 2 = 2 ^7 l7^7(mt) ~ 7i/7(cxt)| 2 (?) 

1 /2 

calculated at the j/-th resonance energy. Another important point to underscore is that the signs of T u 7ext and Mjjc 

1/2 

relative to F ' , int x are synchronized. In all the fits presented below we use positive sign because it gave a better fit. 
Finally the background resonance amplitude can be written as 

r l/2 p l/2 

where Er (bg) is the resonance energy of the background resonance, 

T \{ 2 BG) = \/ 2P lA k P^phl(BG), (9) 



?%G) = V2PiAkc,r c ) lc{BG) , (10) 

and 7 7 (sg) i s a complex radiative width amplitude for the decay of the background resonance to a bound state and 
lc(BG) is the reduced width amplitude of the background resonance for the channel c The total resonance width of 
the background resonance is given by ^(bg) — ^p(BG) + ^u(bg)i where ^ c (bg) is the partial resonance width in the 
channel c. Note that all the energies are in MeV but the astrophysical factor S(E) is in keVb. The dimension of the 
reduced particle width amplitude 7„ c is MeV 1 / 2 but 7„ 7 has dimension MeV^fm 3 / 2 for L = 1. 



IV. ANALYSIS 



Altogether we performed 4 different sets of the fits, each set consists of two fits with the boundary conditions 
B c = S C (E 2 ) and B C {E\). First we present two fits without a background pole. Then we perform 2 unconstrained fits 
called fits A(113), and 2 constrained fits called fits B(113). In these fits we use all 113 data points of [Sj. In addition, 
we performed also two fits A(70) with and without the background resonance using only 70 low-energy data points 
of [B] in the region of the first resonance. For all the fits the channel radii in the proton and a channels, r p = 5.03 
fm and r a — 7.0 fm, and the ANC C = 14.154 fm -1 / 2 are kept the same as in pQ and [TU] and the background 
resonance energy is fixed at Er bg = 5.07 MeV. All other parameters are varied to get a best fit. For the case under 
consideration l p = 0, // = 1, J/ = 0, Jr = 1 and 1 = 1. 

In Fig [l] we demonstrate two best fits to 113 data points of without any background pole with the boundary 
conditions in the channel c = p, a B c = S C (E 2 ) and B c = S c (Ei). The best fit with B c = S C (E 2 ) is achieved at 
E 2 = 0.956 MeV and E 1 = 0.2872 MeV and results in 5(0) = 34.2 keVb with x 2 = 2.6 and total \ 2 = 273.1. For 
the fit with B c = S c (Ei) at Eh = 0.30872 MeV and E 2 = 1.0794 MeV we get 5(0) = 34.6 keVb and x 2 = 2.5 and 
X 2 = 259.9. Parameters of the fits are given in Tables [l| and \H\ As we can see both fits are quite good except for the 
bottom between two resonances and the high energy tail. From Fig. 9(b) [B] we can conclude that to get down % 2 to 
the minimum one really needs to increase significantly the ANC as it has been done in |BJ. 
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FIG. 1: (Color online) The astrophysical S(E) factor for the 15 N(p, 7) ls O reaction. The black squares are experimental data 
from Ref. [6]. The red solid line is our unconstrained fit with the boundary condition B c = S C (E2 = 0.956 MeV), which takes 
into account three interfering amplitudes: two 1~ resonances and non-resonant term. No background pole is included. The 
blue solid line is a similar fit with the boundary condition B c = S c (-Bi = 0.30872 MeV). The ANC is C = 14.154 fm~ 1/2 as in 
all other fits. 



Instead of varying the ANC to decrease \ 2 we perform two unconstrained fits A(113) to the Notre Dame-LUNA 
data [5] adding a background resonance. All the parameters except for the ANC, channel radii and the background 
resonance energy are allowed to vary. First we have searched for the best fit for the boundary condition B c = S C (E2), 
E 2 is the energy of the second level which is taken close to the second resonance energy Er 2 = 0.9594 MeV adopted 
in [5] while the first level is varied to get the best fit. We find the best fit at E 2 = 0.956 MeV and E 1 = 0.1662 MeV. 
From this fit, using the Barker's transformation HUH2], we determine the i?-matrix formal reduced widths for the 
second resonance at resonance energy Er 2 = 0.9594 MeV. After that we can find the observable partial resonance 
widths for the channel c using the standard i?-matrix equation [3] 

VC ~ 1 , 2 dS c , I ' \ 1X ) 

L + <L lvc'^E\E=E Rv 
c f —p,a 

where Er u is the resonance energy of the level v. In the second unconstrained fit A(113), we have searched for the 
best fit with the boundary condition B c = S C (E±), where the energy of the first level E\ is near the first resonance at 
Er x = 0.3104 MeV adopted in [6] while the second level is varied to get the best fit. For this boundary condition we 
find the best fit at Ei — 0.30872 MeV and E 2 = 1.0576 MeV. The formal reduced widths and observable resonance 
widths for the first resonance are determined from this fit by shifting the boundary condition from E\ — 0.30872 MeV 
to the first resonance location at Er x — 0.3104 MeV and using Eq. (11 1. In Fig. [2] we demonstrate the astrophysical 



factors S(E) obtained from these two unconstrained A(113) with fixed ANC C = 14.154 fm -1 / 2 and the background 
resonance included. The red solid line is the fit corresponding to the boundary condition at E 2 = 0.956 MeV with the 
normalized x 2 = 1.84. This fit is practically identical to the one in [B] resulting in S(0) — 39.0 keVb in agreement with 
[S]. For the fit A(113) with the boundary condition at E\ = 0.30872, the blue dotted-dashed line, we obtain x 2 = 1.76 
with S(0) = 37.2 keVb. This fit goes slightly lower than the red line at low energies better reproducing the low-energy 
trend of the data. The magenta solid line represents the non-resonant S(E) factor for the ANC C = 14.154 fm -1 / 2 
which has been used for both fits. This ANC is within the experimental interval 13.86 ± 0.91 fm -1 / 2 determined from 
the 15 N( 3 He, c?) 16 reaction lj. Thus adopting a physical ANC we correctly fix the normalization of the external 
direct capture amplitude and the channel radiative width amplitude, and adding the background pole rather than 
varying the ANC way beyond experimental limits [B] we are able to get the same fit as in |BJ. The obtained formal 
reduced widths and observable resonance widths for the fits A(113) along with the corresponding parameters from [B] 
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FIG. 2: (Color online) The astrophysical S(E) factor for the 15 N(p, 7) ls O reaction. The black squares are experimental data 
from Ref. [B]. The red solid line is our unconstrained i?-matrix fit A(113) with the boundary condition B c — S l c (£ , 2), which 
takes into account four interfering amplitudes: two 1~ resonances, non-resonant term and background resonance at 5.07 MeV. 
The blue dotted-dashed line is our unconstrained i?-matrix fit A(I13) similar to the previous one but with the boundary 
condition B c — S c (Ei). The magenta solid line is the non-resonant S(E) factor for the square of the ANC C = 14.154 fm -1 / 2 . 

and [23] are given in Tables and |n] 

In Fig [3] in the logarithmic scale for both axes, we show the band between upper and low limits of the astrophysical 
factor S(E) obtained from the unconstrained fit A(113) with the boundary condition B c — 5 C (£?2) at E^ ~ 0.956 
MeV for the energy region E < 1.7 MeV. The upper (lower) limit with 5(0) = 40.1 keVb (5(0) = 37.9 keVb) and 
X 2 = 3.0 (x 2 = 2.6) of the band corresponds to the fitting of the experimental data which deviate by 1 a up (down) 
from the center which corresponds to 5(0) = 39.0 keVb with x 2 — 1-84. Thus, taking into account the experimental 
uncertainties given in [BJ, we can conclude that our unconstrained fit A(113) with the boundary condition B c = 5 c (i?2) 
results in 5(0) = 39.0 ± 1.1 keVb. However, the logarithmic scales for both axes show the problem with the fitting at 
low energies, where the fit A(113) with the boundary condition B c — 5 c (i?2) deviates from the experimental trend. 
A similar trend is present in the fit of [BJ. The reason for this trend is that our fits and fit in [BJ have been performed 
minimizing the weighted % 2 with the weights A^ 2 , where Aj is the experimental uncertainty at point i. Since the 
relative experimental uncertainties at low energies are larger than in the region between the two resonances and at 
higher energies, the weighted fit underestimates the importance of the low energy region, which is the most crucial 
for determination of the 5(0) astrophysical factor. 

We can compare the results of the fits A(113) with [BJ and |3]. Note that from different fits presented in [3] we 
choose the fit HH(c), Tables II and III, to the data of [3], because they are pretty close to the data of [BJ in the region 
of the first resonance and the selected fit from [3] resulted in 5(0) = 35.2 keVb, which agrees with our results and 
close to [5] 5(0) = 39.6 ± 2.6 keVb. 

As it has been mentioned in 0], there is significant uncertainty in the values of the proton and a partial widths 
for the second resonance what can be also concluded from compilation 24J. That is why there is no recommended 
values for these widths in [J] after analysis of the 15 N(p, a) 12 C data. One of the problems is that the 15 N(p, a) 12 C 
and 15 N(p, 7) 16 reactions put a limitation on the ratio 7 2 Q ,/7 2 a and the El strength ratio of the second and first 
resonances, which should be equal due to the isospin mixture of two 1 _ resonances. This isospin mixture can be 
written as [2H3| 

fa = a\T = > +P\T = 1 >, 

V>2 = /S|T = 0>-a|T = l>, (12) 
where a 2 + (3 2 — 1. Since the a-particle decays of these resonance in 16 O to the ground state are allowed only due to 
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TABLE I: Resonance parameters. Parameters of the i?-matrix fits to the 15 N(p,7) 16 capture 113 data points [6] along with 
the fitting parameters from [S] and |23]. Fits A(113) and B(113) are our fits. The fitting parameters for the second resonance 
are determined from adopting the boundary condition B c — Sc^E?) at the second energy level E2 = 0.956 MeV, which is near 
the second resonance Ea 2 = 0.9594 MeV, with the first energy level E\ = 0.1662 MeV found from the fit. The channel radii, 
r v — 5.03 fm and r a = 7.0 fm, and ANC C = 14.154 fm -1 / 2 have been used in all the fittings. In all the fittings where the 
background resonance is included the background resonance energy is 5.07 MeV, the proton reduced width amplitude of the 
background resonance is 7 p (bg) — —0.3 MeV 1//2 and the a reduced width amplitude 7 a (sG) = 0.07 MeV 1//2 . In the fit A(113) 
the search for the best fitting has been performed using unconstrained variation of other parameters. Using the Barker's 
transformation [121 122] the fitting parameters are transformed to the ones corresponding to the boundary condition at the 
second resonance energy Er 2 — 0.9594 MeV adopted in [B]. These parameters are given as the fitting parameters for the second 
resonance. The radiative width for the background pole is found to be r 7 ( S G) = 354.9 eV. To determine the fitting parameters 
for the first resonance the boundary condition B c = S c (Ei) has been adopted, where the first energy level E 1 = 0.30872 
MeV. The second energy level found from the fit is E2 = 1.0576 MeV. The fitting parameters are transformed to the ones 
corresponding to the boundary condition at the first resonance E^ = 0.3104 MeV adopted in |6j and are shown in the table 
as the fitting parameters for the first resonance. The radiative width for the background pole is found to be F 7 ( S g) = 360.5 
eV. In the constrained fits B(113) the procedure is the same as described before but two more parameters are fixed. When 
searching for the fitting parameters for the z/-th resonance the proton and a reduced width amplitudes of the resonance v are 
fixed. Their values are taken from the 15 N(p, q) 12 C data fitting parameters [UEI]. in the nrs t constrained fit B(113), which is 
used to determine the parameters for the second resonance, the boundary condition B c = S C (E2) is taken at E2 = 0.956 MeV 
with the first energy level found from the fit E\ — 0.1697 MeV. The radiative width for the background pole is found to be 
r 7 (sG) = 129.3 eV. To determine the parameters for the first resonance we use the boundary condition B c — S C {E\) near the 
first resonance E\ — 0.30872 MeV and found from the fit the second energy level E2 = 1.0576 MeV. The rest is the same as in 
fit A(113). The radiative width for the background resonance is F 7 ( S g) = 283.1 eV. F„ c is the observable resonance partial 
width in the channel c and r„ 7 is the v-th resonance radiative width. 



Reference 


7i% [keV] 


tip [keV] 


7i% [keV] fi 


a [keV] 


Ti 7 [eV] 


72% [keV] f 2p [keV] 


72% [keV] f 2 


ol [keV] T 27 [eV] 


m 


52.8 


0.20 


13.5 


112.0 


33.8 


309.1 110.6 


5.0 


40.6 38.7 


Table II, HH(c) 


355.2 


1.0 


10.6 


85.8 




265.2 98 


5.4 


40.3 


present work, unconstrained fits 
A(113) 


358.8 


1.3 


14.4 


99.4 


7.5 


221.6 82.8 


7.5 


63.2 63.6 


present work, constrained fits 
B(113) 


259.8 


1.0 


13.6 


99.7 


9.3 


268.9 98.1 


6.0 


49.4 54.4 


present work, unconstrained fits 
without background resonance 


353.3 


1.3 


14.1 


98.3 


7.5 


231.4 86.0 


6.9 


58.2 57.6 



the T — components, we have 

- — (13) 

Correspondingly, the strength of the El decays of these resonances to the ground state is entirely determined by the 
T = 1, i.e. 

^ = - (14) 

From the fit in [5] one gets the ratio Jia/j^a = 2.7 and 7i Q /72 Q = 1-96 from [J]. To get the ratio of the El 
intensities we remind that 7 l/7 (i„t) ~< ^\0\^ v ^ int ) > | r <r p , where ip is the ( 15 Np) bound state wave function, O 
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TABLE II: Internal and external radiative width amplitudes for the first and second resonances. The amplitudes for the second 
resonance at Er 2 — 0.9594 MeV (first resonance at E^ — 0.3104 MeV) are determined from the unconstrained fit A(113) and 
constrained fit B(113) for the boundary condition at the second (first) resonance. 



Fits 7l7(int) [MeV 1 / 2 mi 3 - 72 ] 7i 7(Mt) [MeV 1 / 2 fm 3 / 2 ] j^nt) [MeV 1 / 2 fm 3 / 2 ] ~/ 2l(ext) [MeV 1/2 fm 3 / 2 ] 



a 

12J 


n 99 




n in 
u.iy 




in 


0.085 




0.24 




Unconstrained fits A(113) 


0.062 


0.061 + i 0.000059 


0.28 


0.053 + i 0.0041 


Constrained fits B(113) 


0.085 


0.052 + i 0.000054 


0.25 


0.059 + i 0.0044 


Unconstrained fits 
without background resonance 


0.062 


0.060 + i 0.000059 


0.26 


0.055 + i0.0042 




: i i i i i i i i I : 

0.1 1 

E(MeV) 

FIG. 3: (Color online) The astrophysical S(E) factor for the 15 N(p, 7) ie O. The band for the astrophysical factor S(E) obtained 
from the unconstrained fit A(113) with the boundary condition B c — S C (E2). The upper and low limits of the band correspond 
to the fitting of the experimental data, which deviate by 1 a up and down from the center, correspondingly. Note that the 
borders of the band have practically the same particle reduced width amplitudes 7„ c for the first and second levels. The proton 
partial width for the second resonance within the band is T2 P — 82.8 ± 0.6 keV, the a-particle partial width for the second 
resonance F2 Q = 63.2 ± 0.9 keV, the radiative width of the second resonance 1*2 7 = 63.6 ± 2.4 eV and the radiative width of 
the background pole Tj^bg) = 354.9 ± 23.6 eV. The black squares are experimental data from Ref. [6]. 
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is the electromagnetic operator and ip u a n t) is the internal resonant wave function of the v-th resonance given by a 
standing wave satisfying Eq. (12). Then the ratio of the El intensities can be estimated from the ratio of J^, 1 u nt \ 
assuming the dominance of the internal contribution to the electromagnetic transition matrix element. From [6] we get 
72 7 (in*)/Ti7(mt) = °- 66 and 72 7 (mt)/7? 7 ( int ) = 7 - 97 Ir0m SI- If we use for tne E1 intensity ratio of the total radiative 
widths we get from [B] E\ijE\\ — ^-y/fcf 7 )/(ri 7 /fci 7 ) = 0.98, where fc„ 7 = (Er v +e)/hc is the momentum of the 
emitted photon for transition from the resonance v to the ground state with the proton binding energy e. Thus the 
ratio of the a reduced widths from [B| is pretty consistent with findings in [31 [TU] but deviates from [23] , while the 
ratio of the radiative resonance widths is too small compared to all the previous estimations due to too high radiative 
width of the first resonance which was estimated to be around 10 eV [3J[TT1[23]. The proton partial width T2 P = 110 
keV is higher than the previous estimations j|4 s 23 2 1 . Note that we do not include estimations from the analysis of 
the data of [2J. 

The partial widths for the first resonance are better known than for the second one. According to [3] and [23J 
Tip = 1.1 keV and 1.0 keV, correspondingly, and f\ a = 92 ±8 keV and different previous estimations are pretty close 
to these values [24 . Note that all the widths are in the center-of-mass system. That is why too low value of Ti p = 0.2 
keV obtained in [6J is difficult to explain. Our unconstrained fits A(113) is not satisfactory although it better agrees 
with the previous estimations for the first resonance than the fit of [BJ. The quoted value Ti 7 = 12 ± 2 eV in [24] was 
taken from [2J, while [3J obtained from a single level analysis (only the first resonance was included) T l7 = 8 eV, and 
from the two-level analysis Ti 7 = 12.8 eV. Our Ti 7 = 7.5 eV is significantly lower than the corresponding value in [BJ 
and closer to [3JQI]. Our T 27 = 63.6 eV is higher than r 27 = 32 ±5 eV Jl], 38.7 eV [6] and 44 ±8 eV obtained from 
the branching ratio [23] but lower than 88 eV [3J. Our r\ Q = 99.4 keV is in a perfect aggreement with estimation 
92 ± 8 keV [4] and with other estimations [24] while the value Ti Q = 112.0 keV [6] looks beyond of the boundaries of 
the existing estimations. However, our jf a 

h\ a = 1-92 is much smaller than (r 27 /fcf 7 )/(ri 7 /fc? 7 ) = 7.2. 
Due to the above mentioned problems with the unconstrained fits A(113) we performed two constrained fits B(113) 
keeping in mind that with a slightly larger x 2 than for the unconstrained fits we can get more reasonable fitting 
parameters. Once again we did two different fits corresponding to two boundary conditions with parameters given 
in Tables |T| and |H] In these fits, in addition to the fixed the channel radii in the proton and a channels, r p = 5.03 
fm and r a = 7.0 fm, the ANC C = 14.154 fm -1 / 2 and the background resonance energy Er bg — 5.07 MeV, we also 
fix 7„ c , c = p, a, when the boundary condition is chosen near the resonance energy Er v . These reduced widths 
are taken from the analysis of the direct (p,a) data [4] [TO] [IT] and indirect data [10]. First we adopt the boundary 
condition near the second resonance at E 2 = 0.956 MeV with the energy of the first level E\ = 0.170 MeV. For the 
best fit we get \ 2 = 1.93 and 5(0) = 38.8 keVb. After that we set up the boundary condition at E\ — 0.30872 
MeV. For the best fit for the energy of the second level E 2 = 1.0573 MeV we get x 2 = 1.74 [25] with 5(0) = 37.2 
keVb. The parameters given in Tables [IJ and [ITJ are obtained for the boundary conditions at the resonance energies 
adopted in [BJ. Barker's transformation to get the fitting parameters at the energy of the first resonance practically 
didn't change them because of the proximity of our adopted first level E\ = 0.30872 MeV and the first resonance 
energy E^ — 0.3104 MeV adopted in [BJ. In Fig. [4] the S(E) factors are shown for both constrained fits B(113): 
the solid red line represents the fit with the boundary condition at E% = 0.956 MeV and the blue dotted-dashed 
line, which better reproduces the low-energy experimental trend, corresponds to the fit with the boundary condition 
at Ei = 0.30872 MeV. The constrained fits B(113) have parameters which better agree with previous estimations 
[24j than unconstrained fits A(113). In particular, r2 a is lower and T\ p better agrees with previous estimations 
[24] than the width obtained in [6]. The same is true for the observable partial widths for the first resonance which 
agree with previous estimations and the radiative width for the first resonance is in better agreement with |24j . Our 
r2 7 = 54.4 eV, although lower than in the fit A(113), is still high and remains the only problem to get consistency 
with previous estimations [31 [JU1 [23] [23] . Because of that the ratio (r2 7 /fcf 7 )/(ri 7 /fc f 7 ) = 5.0 is still too high 
compared to J 2 a /j 2 a = 2-3 which is fixed in agreement with 15 N(p, a) 12 C (4j [TO] . In Table III we present \ 2 and 5(0) 
astrophysical factors for all our fits and fit from [6]. The uncertainties of our 5(0) factors are obtained by fitting to the 
upper (lower) border of the data obtained by adding (subtracting) the experimental uncertainties to the experimental 
astrophysical factor at each point. A similar procedure has been used to determine the band shown in Fig. [3] As we 
can sec from Table [III] five different fits result in quite stable 5(0) factors ranging in the interval 33.1 < 5(0) < 40.1 
keVb. The unconstrained fit A(113) and constrained fit B(113) with B c = S c (Ei) give the minimum \ 2 among all 
our fits with 5(0) = 37.2 ± 1.0 keVb. which overlaps with the result reported in [B]. However, as we have discussed, 
the constrained fits B(113) yield fitting parameters including a correct ANC value, which are more consistent with 
the previous estimations. Besides, our both best fits better reproduce the low-energy slope of the S{E) astrophysical 
factor than the unconstrained fit A(113) with B c — S C (E 2 ), see Fig [3j and the fit of [BJ, which trend away from 
the low energy data, and a lower value of the 5(0) is quite plausible when lower energy data will be available. We 
included also into the list of the fits two fits performed without a background resonance because their \ 2 deviate from 
the minimum % 2 by < 1. These fits provide the lowest 5(0) better reproducing the low energy behavior of the S(E) 



i — 1 1 1 — i — • — i 1 — i — • — i — • — r 




1 r 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 
E(MeV) 



FIG. 4: (Color online) The astrophysical S(E) factor for the 15 N(p,7) 16 C\ The red solid line is the constrained fit B(113) with 
the boundary condition B c — S C (E2 = 0.956 MeV), the blue dotted-dashed line is the constrained fit B(113) with the boundary 
condition B c = S c (Ei = 0.30872 MeV). The black squares are experimental data from Ref. [6]. The ANC is C = 14.154 
fm- 1 / 2 . 



TABLE III: x 2 and S(0) astrophysical factors for the 15 N(p, 7) 16 capture process obtained from our fits and from the fit in 



Fits 


x 2 


S(0) [keVb] 


Ref. [6] 


1.80 


39.6 ±2.6 


Fit A(113), B c = S C {E 2 ) 


1.84 


39.0 ± 1.1 


Fit A(113),S C = 5 c (£i) 


1.76 


37.2 ± 1.0 


Fit B(113),B C = S C (E 2 ) 


1.93 


38.8 ± 1.1 


Fit B(113), B c = Sc(Ei) 


1.74 


37.2 ± 1.0 


Fit without background resonance, 






B c = S C (E2) 


2.58 


34.1 ± 1.0 


Fit without background resonance, 






B c = S c (Ei) 


2.45 


34.6 ± 1.0 



10 



400 




I ' 

0.1 0.15 0.2 0.25 0.3 0.35 0.4 

E(MeV) 



FIG. 5: (Color online) The astrophysical S(E) factor for the 15 N(p, 7) ie O reaction in the low-energy region including the first 
resonance (70 data points). The black squares are the experimental data from Ref. [6]. The red solid line is the unconstrained 
fit A(70) for the boundary condition B c = S c (Si = 0.30872 MeV); for this fit E 2 = 1.056 MeV, y 2 -f(int) = 0.056 MeVfm 3/2 , 
J2~,(ext) = 0.0607 + i 0.000059 MeVfm 3/ ^ 2 . The blue dotted-dashed line is the constrained fit B(70) with the boundary condition 
B c = S c (Ei = 0.30872MeV); for this fit E 2 = 1.05 MeV, -y 2l{int ) = 0.082 MeVfm 3/2 , -y 2l{£Xt) = 0.052 + i 0.000051 MeVfm 3/2 . 
The ANC is C = 14.154 im~ 1/2 . 



TABLE IV: Resonance parameters. Parameters of the unconstrained A(70) and constrained B(70) fits to the 15 N(p, 7) ls O 
70 data points [6]. Physical meaning of the parameters and the procedure are similar to the one described in the caption for 
Table [I] Note that for both fits A(70) and B(70) we present the parameters only for the first resonance because the boundary 
condition is adopted at the energy of the first level. 





Fit 


7l % [keV] fi, [keV] 


7 L [keV] fia [keV] Fi 7 [eV] 




fit A(70), 






B c 


= S c (Ei = 0.30872 MeV) 


358.9 1.4 


14.3 98.9 6.9 




fit B(70), 






Be 


= S c (£i = 0.30872 MeV) 


260.2 1.0 


13.6 99.7 9.0 



factor than the ones with higher 5(0) . 

To increase the weight of the low-energy points we present also two fits to 70 low-energy data points in the region 
of the first resonance rather than to all 113 data points: the unconstrained fit A(70) and the unconstrained fit B(70) 
with the boundary condition B c = S c (Ei). New fits to 70 data points are shown in Fig. [5] and parameters are given 
m Tables [!V| and |y} We use the the logarithmic scale for the energy axis and linear scale for the S(E) factor to see 
more clearly the low-energy behavior of the astrophysical factors. 

It is worth mentioning that both fits to 70 data points favor lower value of the 5(0) factors than in i.e. the 
same tendency which we have observed for the corresponding fits to 113 data points. The parameters of both fits are 
similar to the ones of the fits for 113 data points. 
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TABLE V: X 2 an d S(0) astrophysical factors for the 15 N(p, 7) ie O capture process obtained from the unconstrained fit A(70) 
and constrained fit B(70). 





Fits 




x 2 


5(0) [keVb] 


Fit A(70), B c 


= S c (£i = 


: 0.30872MeV) 


1.50 


37.7 


Fit B(70), B c 


= S e (Ek = 


0.30872MeV) 


1.51 


37.1 




V. 


SUMMARY 







Determination of the 5(0) factor for the 15 N(p, 7) 16 radiative capture is one of the goal of our fits. Although 
new measurements of this reaction [S] is a real success and a very important contribution to study of this reaction, 
we believe that it would be difficult to give a more accurate 5(0) value than the range 33.1 < 5(0) < 40.1 keVb 
determined from our fits without further measurements down to lower energies than those achieved in [6]. To get 
more accurate uncertainties of the 5(0) factor a better estimate of energy uncertainties would be also useful. From our 
fits we determine the interval of the astrophysical factors at the effective energy E = 23.44 keV 36.0 keVb < S(E = 
23.44 keV) < 44.46 keVb. Assuming the astrophysical factor 84.1±5.9 MeVb for the competing reaction 15 N(p, cv) 12 C 
[TU] we find that for every 2084+g2g cycles of the main CN cycle one CN catalyst is lost due to the 15 N(p, 7) 16 
reaction. 

But what is even more important is the question whether the minimum of x 2 is always an acceptable fit. Definitely, 
the answer is yes if our knowledge about the reaction model is complete. But it is assumed that the best fit is 
achieved under constrained variations of the fitting parameters within the accepted boundaries obtained from the 
available physical information. This question elevates when the input physics is not complete. It is another goal of 
our analysis to demonstrate that, due to our incomplete knowledge of the reaction model, it is not canonical that a 
fit, which provides minimum of % 2 , is the best from the point of view of physics. We have demonstrated here that it is 
possible to achieve the same or even better fit and similar final 5(0) factors as in [B] by adopting the ANC measured 
from the transfer reaction rather then using an unconstrained variation of the ANC. But, even if the ANC is fixed 
within the experimental boundaries, the question remains about the interpretation of other fitting parameters. We 
have demonstrated problems with the interpretation of the parameters of the fits A and fit in [B]. Trying to improve 
interpretation we fixed some parameters, which are available from the analysis of the 15 N(p, a) 12 C reaction, and 
we are able to achieve even better fit than in [5J and better agreement of the fitting parameters with the previous 
measurements of the 15 N(p, 7) 16 0, 15 N(p, a) 12 C and 15 N(p, p) 15 N processes [24]. However, still even our constrained 
fits are not fully satisfactory because we got too high value of the radiative width of the second resonance if we assume 
that the radiative width of the first resonance is ~ 10 eV. This issue remains to be resolved. 
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